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ABSTRACT 

PLOT3D is a package of FORTRAN subprograms to draw three-dimensional surfaces 
of the form z = f(x,y). The function f and the bounding values for x and y are the 
input to PLOT3D. The surface thus defined by PLOT3D may be drawn after arbitrary 
rotations. Output is by off-line incremental plotter or on-line microfilm recorder. 
PLOT3D is completely described along with its limitations, possible future modifications, 
and potential uses. Sample output and listings of the subprograms are included. Written 
entirely in FORTRAN IV, PLOT3D is readily adaptable to other hardware. 

STAR Category 08 


ii 


views are not perspective. PLOT3D is designed primarily to portray functions of two 
independent variables. The functions must be piecewise continuous and bounded. 

This report contains a complete description of the PLOT3D package. In addition, 
potential applications are presented; examples of PLOT3D output are included; and the 
limitations of PLOT3D are discussed. 

PLOT3D is implemented for both the IBM 7044/7094 Direct Couple System - 
California Computer Products Digital Incremental Plotter, controller model 780, printer 
model 765; and the IBM 360/67 - Control Data Corporation microfilm system model 280 
(referred to at Lewis as the DD80). The subprograms for both are explained, and there 
is a sample user’s program. Finally, possible future modifications are mentioned. 


POSSIBLE APPLICATIONS FOR PL0T3D 

It can be extremely difficult to get a ’’feel” for the behavior of complicated functions 
or the shape of empirical surfaces from tabulated values or conventional contour plots. 
Since plotting such surfaces is PLOT3D's purpose, PLOT3D may be of use anywhere 
these surfaces arise. 

In engineering and physics, functions of two variables arise often, both experimen- 
tally and theoretically. PLOT3D can portray their behavior. For example, the phase 
diagrams and equations of state of thermodynamics can be studied. In mechanics, for 
example, energy surfaces and potential wells can be drawn (see fig. 1). Data points 
from experiments can be curve-fitted and drawn in three dimensions (see fig. 2, also 
ref. 8). 

Partial differential equations in two independent variables occur frequently - for ex- 
ample, in heat flow, electricity and magnetism, fluid dynamics, and quantum mechanics. 
Maps of their analytic and numerical solutions are easy to make with PLOT3D (see 
fig. 3). 

In numerical analysis and analytic geometry, problems arise which can be solved 
with the aid of PLOT3D: The behavior of functions near singularities can be pictured 
(see figs. 4(a), (b), and (c)); approximate roots of equations can be located; and the form 
of surfaces after rotations, contraction, and deformations can be pictured. Numerical 
solutions to equations of any sort can be curve-fitted and then mapped by PLOT3D. 


LIMITATIONS OF PL0T3D 

PLOT3D is designed to draw only functions; that is, for every point in the x and y 
ranges there must be a unique z value. It cannot draw, for example, a sphere. 
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PL0T3D - A PACKAGE OF FORTRAN SUBPROGRAMS TO DRAW 
THREE-DIMENSIONAL SURFACES 
byR. Bruce Canright, Jr., and Paul Swigert 
Lewis Research Center 

SUMMARY 

PLOT3D is a package of FORTRAN subprograms to draw three-dimensional surfaces 
of the form z = f (x, y) . The function f and the bounding values for x and y are the 
input to PLOT3D. The surface thus defined by PLOT3D may be drawn after arbitrary 
rotations. Output is by off-line incremental plotter or on-line microfilm recorder. 

Plot3D is completely described along with its limitations, possible future modifica- 
tions, and potential uses. Sample output and listings of the subprograms are included. 
Since PLOT3D is written entirely in FORTRAN IV, it is readily adaptable to other hard- 
ware. 


INTRODUCTION 

Several packages for plotting three-dimensional surfaces have been developed else- 
where. Each package has its strengths and limitations. Some of these programs are in- 
tended primarily to show the spatial relations among solids. This problem was perhaps 
studied first by Roberts, who considered convex polyhedra (ref. 1). Appel has developed 
powerful methods for drawing solids bounded by planes (ref. 2). Weiss has dealt with 
surfaces which are plane or quadric (ref. 3). 

Other packages draw surfaces made up of data points. Perspective drawings can be 
made from a given viewpoint (refs. 4 and 5). Branstetter has designed a program to 
draw contour maps (ref. 6), and there are similar programs available commercially. 
Nelson has developed routines which, within certain restrictions, will plot both arrays of 
data and functions of two variables (ref. 7) . All packages are machine dependent to some 
extent. 

PLOT3D will plot any function of the form z = f(x, y). It uses a given recipe to gen- 
erate views of the surface after arbitrary rotations about the three spatial axes. The 



The functions must be in rectangular coordinates. Also, the function must be ex- 
pressed explicitly in the form z = f(x, y). Where such a form is inconvenient or impos- 
sible, this restriction can sometimes be overcome by curve fitting methods. 

Every line lying on the surface is drawn, not just those which are "visible" or "in 
front. " In a sense, the views are projections of wire frame models onto a plane. They 
are not like wood models. 

The output views may be difficult to study quantitatively. It requires careful choice 
of the input parameters to read values off the surface plots. 


DESCRIPTION OF PL0T3D PACKAGE 

The method used by the routines is simple. The given x and y ranges are divided 
to form a uniform rectangular grid. The values of the supplied function at the grid points 
(x, y) are calculated and stored; this defines the surface. 

The surface is portrayed by connecting successive (y, z) points with straight-line 
segments for each x value on the grid and, in turn, connecting successive (x, z) points 
for each fixed y value on the grid. These lines are then projected by parallel projection 
onto the fixed yz -plane for plotting. 

The function is evaluated only once, before the first view. Views of the surface are 
made after rotating it about the fixed spatial axes - first around the fixed x-axis through 
angle a (see fig. 5(a)), then around y through /3 (fig. 5(b)), and finally around z 
through y(fig. 5(c), see also ref. 10). For this order it is seen that the matrix which 
performs the desired rotation is simply the product CBA where A, B, and C are the mat- 
rices defined in figure 5. The grid points and their function values are transformed and 
redrawn. Note that the three axes of rotation are fixed and that the order of rotations is 
crucial, because matrix multiplication is not in general commutative: for example, the 
rotation ABC does not give the same view as CBA. 

Input parameters allow considerable flexibility in drawing the surface. The x and 
y ranges are determined by supplying maximum and minimum x and y values. The 
number of grid lines in each direction is supplied, as well as the number of points on 
each grid line. The more grid lines there are, the finer the mesh is. This fine mesh 
allows shading of areas with small surface gradients. Because straight-line segments 
are drawn between the points, the more points that are used, the better the representa- 
tion of curvature is. The number of points is limited by the computer storage space and 
by the computation time. Finally, the surface may be rotated about the three fixed axes, 
as discussed earlier. 

A general flow chart of the subprograms which make the drawings is given in fig- 
ure 6. As the left side of this figure shows, SUBROUTINE PLOT3D sets up the rectan- 
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gular grid and uses the supplied function to get the z values. At this point the figure 
axes coincide with the fixed axes; that is, the figure yz plane coincides with the fixed yz 
plotting plane. 

ENTRY ROTATE to PLOT3D calls subroutines which draw the surface after arbi- 
trary rotations about the three fixed axes, as the right side of figure 6 shows. A de- 
scription of the routines called by PLOT3D and ROTATE follows. 

SUBROUTINE SCALE, called once for each surface to be drawn, locates the mini- 
mum and maximum z values. SCALE then scales x, y, and z to fit into a cube 3. 5 
plotting units on a side, with the plotting plane defined as a square with 10 plotting units 
on a side. 

The variables x, y, and z may be scaled in one of two ways. Each variable may 
be scaled according to the maximum magnitude it takes on, or all three variables may be 
scaled according to the largest of these three maxima. The first way of scaling fits the 
surface into a cube and may distort it considerably. The second way of scaling pre- 
serves the original relative magnitudes, but may shrink one or two axes so small as to 
be invisible. The bounding values are placed in COMMON/MAXES/XMIN, XMAX, YMIN, 
YMAX, ZMIN, ZMAX. 

SUBROUNTINE AXIS, called once for each surface view, sets up and draws the 
three figure axes. These axes help to show the orientation of the surface with respect to 
the fixed spatial axes (see figs. 5 and 7). AXIS calls PHI to rotate these figure axes. 

SUBROUTINE DRAW, called once for each view plotted, draws the line segments 
connecting (y,z) grid points. DRAW, like AXIS and WRITE, uses plotting routines which 
are machine dependent. 

SUBROUTINE WRITE, called once for each view plotted, labels the figure with the 
bounding values from COMMON/MAXES/XMIN, XMAX, YMIN, YMAX, ZMIN, ZMAX. 

This subroutine also labels the figure with the accumulated angles a, /3, and y, and 
with an alphanumeric array supplied through COMMON/LABEL/ as explained in the sec- 
tion USING PLOT3D. The array in COMMON/ LABEL/ is initially filled with blanks in a 
BLOCK DATA subprogram. 

SUBROUTINE TRNMAT sets up the rotation matrix M (M = CBA where the matrices 
A, B, C are defined according to fig. 5) by using the supplied angles of rotation a } 
and y. TRNMAT is called once for every rotated view, that is, every call to ROTATE. 
The rotation matrix is placed in COMMON/MATRIX/TMAT (3, 3). 

SUBROUTINE PHI rotates the surface curves and figure axes by multiplying the po- 
sition vector of each point by the rotation matrix in COMMON/MATRIX/TMAT (3,3). It 
is called three times for each rotated view (see fig. 6). 

The output devices and subprograms are dependent on the computer used to calculate 
and store the points to be plotted. PLOT3D has been written for both the IBM 7044/94 
DCS and the IBM 360/67. On the 7044/94, PLOT3D uses the CalComp programs 
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SYMBOL, NUMBER, and PLT770 (ref. 8). Output is written onto tape which is then used 
by the off-line CalComp plotter. On the 360/67, PLOT3D uses the plotting routine 
LRCURV and several associated routines developed at Lewis for the on-line microfilm 
recorder. Presently, a positive film is produced, which gives white -on-black prints. 

Subroutine listings for the IBM 7044/94 appear in appendix A, while listings of the 
subroutines for the IBM 360/67 are given in appendix B. The functions of the plotting 
programs used have been explained in comments cards so that similar routines at other 
facilities might be substituted. 


USING PL0T3D 

The first call must be CALL PLOT3D(XMIN, XMAX, YMIN, YMAX, S, NXPTS, NYPLS, 

NXPLS, NYPTS, FUNC, CUBE) 

where 

XMIN minimum x value 

XMAX maximum x value 

YMIN minimum y value 

YMAX maximum y value 

S name of the array used by PLOT3D to hold all the coordinates of points to be 

plotted (It must be dimensioned at least 3*(NXPTS*NYPLS+NYPTS*NXPLS) . ) 

NXPTS number of points in each plane y = constant (These are connected by lines to 
represent a curve on the surface.) 

NYPLS number of planes y = constant 

NXPLS number of planes x = constant 

NYPTS number of points lying in each of these planes 

FUNC name of the function z = f(x,y) to be plotted (It must appear in an EXTERNAL 

statement in the user's calling program. The user must supply a function 
subprogram of the form FUNCTION FUNC(x,y).) 

CUBE LOGICAL variable (When CUBE = . TRUE. , the three dimensions are scaled 
independently in order to fit the surface into a cube. When CUBE = . FALSE. , 
the surface is not fit into a cube. See the discussion on SUBROUTINE SCALE, 
p. 4.) 

The preceding parameters, except for S, FUNC, and CUBE are illustrated in figure 7. 

Views of the surface may now be made by the call CALL ROTATE (A, B, C, DRAWME) 
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where 

A angle of rotation about the fixed x axis, deg 

B angle of rotation about the fixed y axis, deg 

C angle of rotation about the fixed z axis, deg 

DRAWME LOGICAL variable set to . TRUE, if the particular view is to be plotted, 

or set to . FALSE, if this view is not to be plotted 

Positive angles are counterclockwise rotations, while negative angles are clockwise ro- 
tations. The surface is not returned to its original position after rotation. Labels may be 
placed below the views by placing an alphanumeric array in the common block/ LA BEL/. 
For the 7044/94, COMMON/LABEL/LAB(12), this array may contain up to 12 words of 
6 characters each; for the 360/67, COMMON/LABEL/LAB(15), the array may contain 
up to 15 words of 4 characters each 

A sample calling program with function subprogram for the 7094/44 follows. This 
program produced figures 4(a), (b), and (c). Note that for use on the 360/67 only the 
array LAB is changed. 


EXTERNAL POLES 
LOGICAL CUBE, DRAWME 
DIMENSION LABAI 12 ) ,S (6000) 

COMMON/LABEL/ LABI 12) 

DATA LABA /2*6H ,6HW = CA, 6HBSI1/Z, 6H + l/I, 6HZ-1)**, 

* 6H2) . Z, 6H = X +, 6H I *Y 3*6H / 

CUBE = .TRUE. 

DO 1 1=1,1? 

1 LABI I ) = LABAI I ) 

CALL PL0T3DI-3.0,3.0,-3.0,3.0,S,51 ,31,7,51 , POLES, CUBE) 

DRAWME = .FALSE. 

CALL ROTATE I 0.0, 0.0, 4 5.0, DRAWME ) 

DRAWME = .TRUE. 

DO 2 1=1,3 

CALL ROTATE I 0.0, 2 0.0, 0.0, DRAWME) 

2 CONTINUE 
STOP 

END 
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FUNCTION POLES ( X . Y ) 

C FUNCTION WITH TWO POLES OF DIFFERENT ORDER 

COMPLEX Z 

IF(Y.NE.O.O) GO TO 1 
I F ( X.NE.0.0 .AND. X.NE.l.O) GO TO 1 
FA = ?0. 

GO TO 2 

1 Z = CMPL X ( X » Y ) 

FA = CABS ( 1 »0/ Z + 1 .0/ ( Z-l .0 ) **2 ) 

I F ( FA.GT.20. ) FA = 20. 

2 POLES = FA 
RETURN 
END 

The PLOT3D package requires about 950 jq storage locations excluding the plot rou- 
tines and the S array of coordinates. Execution time is roughly 1 to 3 seconds per view 
for the examples given in figures 1 to 4 and ~7500 points on the 360/67. There are no 
error returns in PLOT3D. For the 7044/94 computer, the execution time is less, but 
each view requires 1 to 2 minutes of plotter time. 


CONCLUDING REMARKS 

A package of FORTRAN subprogram has been designed to draw three-dimensional 
surfaces by means of either the 7044/7094 DCS and off-line CalComp plotter, or the 
360/67 and on-line CDC microfilm recorder. The package is readily adaptable to other 
hardware. It is intended primarily for functions of two independent variables and is not 
much use in drawing, for example, complex three-dimensional solids. 

The PLOT3D package can be of use in many areas of mathematics, physics, and en- 
gineering. It has already had successful applications here at Lewis. 

PLOT3D could have new uses when new devices are available with the 360/67. With 
a scope and a conversational mode of operation, the user could draw, rotate, and alter 
many surfaces in a short time. Motion pictures have been made, which adds the fourth 
dimension - time. 

Future modifications could increase the powers of PLOT3D. The projection of views 
onto the fixed plane could be made a point projection; that is, perspective views could be 
made. Functions in coordinate systems other than rectangular could be plotted. Parts 
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of lines which are "hidden" could be determined and left off the view. Probably the only 
method to accomplish this would be a point-by-point visibility check, which is a time- 
consuming process. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 27, 1968, 

129-02-01-05-22. 
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APPENDIX A 


Routines for Direct Couple System 


The following are listings of the subroutines contained in the PLOT3D package for 
the IBM 7044/94 DCS. These routines are used with an off-line CalComp digital incre- 
mental plotter. The standard CalComp plotting routines used are explained in the com- 
ments cards. 
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SUBROUTINE PL0T3D 
PURPOSE 

DRAW AND ROTATE THREE DIMENSIONAL SURFACES WITH THE CALCOMP 
DIGITAL PLOTTER 

USAGE 

CALL PLOT3D(XMIN,XMAX,YMIN,YMAX,S,NXPTS,NYPLS,NXPLS,NYPTS, 

Ft CUBE ) 

CALL ROTATE ( ALPHAt BETA, GAMMA tDRAWME) 

COMMON /LABEL/LAB ( 12 ) 

DESCRIPTION OF PARAMETERS 

XMIN -MINIMUM •X' VALUE TO BE PLOTTED 

XMAX -MAXIMUM »X« VALUE TO BE PLOTTED 

YMIN -MINIMUM «Y« VALUE TO BE PLOTTED 

YMAX -MAXIMUM 'Y 1 VALUE TO BE PLOTTED 

S -ARRAY OF LENGTH 3* ( NXPTS*NYPLS+NXPL S*NYPTS ) 

NXPTS -NUMBER OF POINTS DESCRIBING LINES IN PLANES PARALLEL 
TO THE XZ PLANE 

NYPLS -NUMBER OF PLANES PARALLEL TO THE XZ PLANE 

NXPLS -NUMBER OF PLANES PARALLEL TO THE YZ PLANE 

NYPTS -NUMBER OF POINTS DESCRIBING LINES IN PLANES PARALLEL 
TO THE YZ PLANE 

F -FUNCTION SUBPROGRAM SUPPLIED BY THE USER TO DESCRIBE 

THE SURFACE TO BE DRAWN. IT MUST BE OF THE FORM 
FUNCTION F ( X t Y ) • 

CUBE -LOGICAL VARIABLE SET TO .TRUE. IF THE THREE 

DIRECTIONS ARE TO BE SCALED INDEPENDENTLY OF EACH 
OTHER, AND SET TO .FALSE. IF THE THREE DIRECTIONS 
ARE TO BE SCALED DEPENDENTLY 

ALPHA -ANGLE OF ROTATION ABOUT THE X AXIS GIVEN IN DEGREES 

BETA -ANGLE OF ROTATION ABOUT THE Y AXIS GIVEN IN DEGREES 

GAMMA -ANGLE OF ROTATION ABOUT THE Z AXIS GIVEN IN DEGREES 

DRAWME -LOGICAL VARIABLE SET TO .TRUE. IF THE ROTATED 

VIEW IS TO BE DRAWN, AND SET TO .FALSE. IF THIS 
VIEW IS NOT TO BE DRAWN. 

LAB -ARRAY OF BCD CHARACTERS USED TO LABEL THE DRAWING. 

(THIS ARRAY NEED NOT BE SUPPLIED) 

REMARKS 

A CALL TO PLOT3D DEFINES AND SCALES THE SURFACE. 
SUBSEQUENT CALLS TO ROTATE ROTATE IT TO NEW FRAME OF 
REFERENCE. THE SURFACE IS NOT RETURNED TO THE ORIGINAL 
POSITION AFTER A ROTATION. THE SURFACE IS ROTATED ABOUT 
THE STATIONARY XYZ AXES IN THIS ORDER+ FIRST ABOUT THE *X* 
AXIS THEN ABOUT «Y* AND FINALLY ABOUT • Z 1 . 
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SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 
PLT770 (A CALCOMP SUBROUTINE) 

NUMBER (A CALCOMP SUBROUTINE) 

SYMBOL (A CALCOMP SUBROUTINE) 

DRAW 

WRITE 

TRNMAT 

PHI 

SCALE 

AXIS 

NOTE — PLT770 MOVES THE PEN TO (X,Y) ON THE PLOTTING 
SURFACE, WITH THE PEN UP OR DOWN. 

NUMBER DRAWS A NUMBER AT (X,Y). 

SYMBOL DRAWS A BCD SYMBOL AT (X,Y). 

A FUNCTION SUBPROGRAM SUPPLIED BY THE USER TO 
DESCRIBE THE SURFACE TO BE DRAWN. THE NAME OF THIS 
FUNCTION SUBPROGRAM MUST APPEAR IN AN EXTERNAL 
STATEMENT OF THE USER'S CALLING PROGRAM. 


****** ************* ******** ********** ********** ****** * ************ 

SUBROUTINE PLOT3D < XMI N , XMAX , YM I N , YMAX , S , NXPT S ,NYPL S , 

* NXPLS,NYPTS,F,CUBE> 

LOGICAL CUBE 
DIMENSION S ( 1 ) 

DATA IK/O/ 

IF ( IK.NE.O) GO TO 1 
CALL WRITE 
I K= I 

1 FNX 1 = NXPTS-1 
FNY I = NYPLS-1 
FNX2 = NXPLS-1 
FNY2 = NYPTS-1 
LI = NXPTS*NYPLS 
L2 = NXPLS*NYPTS 
N1 = L 1 
N2 =N1 +L I 
N3=N2+L 1 
N4=N3+L2 
N5=N4+L2 
DX=XMAX— XMI N 
DY=YMAX— YM I N 
C 
C 

DO 2 I=1,NXPTS 
DO 2 J=1 ,NYPLS 

I I = NX PTS* ( J-l ) + I 

II = 1 I 
I2=N1 + I I 
I 3=N2+ I I 

S< II )=XMIN+DX*FLOAT< I-D/FNX1 
S ( 1 2 ) =YM I N+DY*FLOAT (J-ll/FNYl 
2 S( I3)=F(S( II ) ,S( 12) ) 

C 

c 

DO 3 I = 1 , NXPL S 
DO 3 J=1,NYPTS 
II = NXPLS*(J-1) + I 
I 1 =N3+ 1 I 
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I 2=N4+ I I 
I 3 = N5+ 1 I 

S( II )=XMIN+DX*FLOAT{ 1-1 )/FNX2 
S( I 2 ) = YM I N+DY*FL 0 AT { J-l )/FNY2 
S(I3)=F(S!I1),S(I2)) 


N1 =N1 + 1 
N2=N2+ 1 
N3=N3+ 1 

N4=N4+ 1 

N5=N5+ 1 

CALL SCALE!XMIN,XMAX,YMIN,YMAX,S(1) ,S(N1 ) , S ( N2 ) , NXPT S , NY PL S , 
* S(N3) ,S(N4) ,S(N5) ,NXPLS»NYPTS,CUBE > 


CALL AXIS(0, .FALSE. ) 


ASUM=0.0 

BSUM=0.0 

CSUM=0.0 

RETURN 


ENTRY ROTATE (ALPHA, BET A, GAMMA , DRAWME ) 

LOGICAL DRAWME 

CALL PLOTIO. 0,0. 0,-3) 

CALL PLOTCO. 0,5. 0,-3) 

CALL TRNMAT ( AL PHA , BET A , GAMMA ) 

CALL PHI (SU),S(N1),S(N2) ,NXPTS,NYPLS ) 
CALL PHI <S!N3),S!N4),S!N5),NXPL$,NYPTS) 
CALL AXIS( 1, DRAWME) 

ASUM=ASUM+ALPHA 


BSUM=BSUM+BETA 

CSUM=CSUM+GAMMA 

I F{ .NOT .DRAWME ) RETURN 

CALL DRAWS ( S ( 1 ) , S ( N1 ) , S ( N2 ) , NXPTS , NYPLS ) 
CALL DRAW(S!N3> ,S(N4) ,S!N5) ,NXPLS,NYPTS) 


CALL WRITES! ASUM, BSUM , CSUM ) 
CALL PLOT! 15. 0,-5. 0,-3) 
RETURN 
END 


BLOCK DATA 

COMMON/LABEL/ LAB! 12 ) 
DATA LAB /12*6H 
END 
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SUBROUTINE DRAW ( X , Y , Z , NX , N Y ) 

DIMENSION X ( NX » NY ) * Y ( NX » NY ) , Z(NX,NY) 
DO 1 1=1, NX 

CALL PLOT (Y(I,1),Z(I,1),3) 

DO 1 J=2 , NY 

CALL PLOT (Y(I,J),Z(I,J),2> 

1 CONTINUE 
RETURN 

ENTRY DRAWS ( X, Y , Z ,NX , N Y ) 

DO 2 J= 1 * NY 

CALL PLOT ( Y( 1, J) ,Z ( 1,J ) ,3) 

DO 2 1=2, NX 

CALL PLOT ( Y( I , J) ,Z ( I , J) ,2 ) 

2 CONTINUE 
RETURN 
END 
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SUBROUTINE WRITE 

COMMON /MAXES/ XMI N , XMAX , YMIN , YMAX , ZM I N , ZMAX 
COMMON /XCPIDX/ NBL , RUNNO,CP I D ( 8 ) 

COMMON /LABEL/ LABI 12) 

DATA BLOCK, ALPHA, BETA, GAMMA/ 5HBLOCK , 6HAL PHA=, 5HBET A= , 6HGAMMA=/ 
DATA HXMIN,HXMAX,HYMIN,HYMAX,HZMIN,HZHAX/5HXMIN=,5HXMAX=,5HYMIN 
1 HYMAX= *5HZMIN=*5HZMAX=/ 

CALL PLOT (O.O, O.O, -3) 

CALL PLOT (0.0,5. 0,-3) 

CP I D ( 5 ) =RUNNO 

CALL SYMBOL ( 0.0, -4.0, .25, CP I D, 90 .0, 30 ) 

CALL SYMBOL ( 0.0 ,-5 . 0,0 .05 , BLOCK ,0. 0 , 6 ) 


NBL1=NBL-1 

PNBL=NBL1 

CALL NUMBER ( 0. 3,-5 .0, 0. 05 , FNBL , 0 .0 ,-l ) 

WRITE (6,2) NBL1 

CALL PL0T(9. 0,-5. 0,-3) 

RETURN 

ENTRY WRITES(A1,B1,C1 ) 

A=AM0D(A1, 360.0) 

B=AM0D(B1, 360.0) 

C=AM0D(C1 ,360.0) 

CALL SYMBOL ( 0.0, -5 .0, 0. 05 , BLOCK , 0. 0, 6 ) 
FNBL=NBL 

CALL NUMBER ( 0. 3 ,-5 .0 , 0.05 , FNBL , 0 .0,-1 ) 
WRITE (6,2) NBL 
IF (LAB(l).EO.ISTR) GO TO 1 
CALL SYMBOL ( -3 .0 ,-4.0 , 0. 1 , LAB , 0.0 , 72 ) 
CALL SYMBOL ( -3 . 0 ,-4.25 ,0 . 1 , ALPHA , 0. 0, 6 » 
CALL NUMBER ( -2 .4,-4. 25 ,0. 1 , A ,0.0, 5 ) 

CALL SYMBOL ( -3 . 0,-4. 5 , 0. 1 , BET A, 0 .0, 6 ) 
CALL NUMBER ( -2 .4,-4. 5,0. 1 ,B, 0.0, 5) 

CALL SYMBOL ( -3 .0,-4. 75 ,0. 1 , GAMMA , 0. 0, 6 ) 
CALL NUMBER (-2 .4,-4. 75 ,0. 1 ,C ,0.0,5 ) 

CALL SYMBOL ( -1 .0,-4. 2 5 , 0 . 1 , HXMI N, 0. 0, 5 ) 
CALL NUMBER ( -0 .4 ,-4.25 ,0. 1 , XM I N ,0 . 0, 5 ) 
CALL SYMBOL ( -1 .0 ,-4.5 ,0 . 1 , H YM I N, 0 .0, 5 ) 
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c 

2 


CALL NUMBER 
CALL SYMBOL 
CALL NUMBER 
CALL SYMBOL 
CALL NUMBER 
CALL SYMBOL 
CALL NUMBER 
CALL SYMBOL 
CALL NUMBER 
RETURN 


FORMAT ( 120X,5HBL0CK,I5) 
END 



SUBROUTINE SCALE ( XM I N , XMAX , YM I N , YMAX , X 1 , Y 1 , Z 1 , NX 1 , NY 1 , X 2 , Y2 , 12 , NX 
12,NY2,CUBE) 

DIMENSION X 1 ( NX 1 , NY 1 ) , Yl(NXl f NYl), Z1(NX1,NY1), X2(NX2,NY2>, Y2(N 
1X2, NY2) t Z2(NX2,NY2> 

COMMON /MAXES/ XM I N1 T XM AX 1 , YM I N1 , YMAX 1 , ZM I N , Z MAX 

LOGICAL CUBE 

REAL MAXX,MAXY ,MAXZ 

XMIN1=XMIN 

XMA X 1 = XMAX 

YMIN1=YMIN 

Y MAX 1 = YMAX 

MAXX=( XMAX-XMIN)/3.5 
MAXY = ( YMAX-YMIN) /3.5 
ZMAX=Z1<1,1) 

Z M I N= Z 1 ( 1 , 1 ) 

DO 1 1=1, NX1 

DO 1 J = 1 , NY 1 

Z MA X = AMA X 1 ( Z1 ( I , J ) , ZMAX ) 

1 ZM I N= AM I N 1 ( Z 1 ( I , J) , Z M I N ) 

DO 2 1=1, NX2 

DO 2 J = 1 , NY2 

Z MA X=AMAX1 (Z2< I , J ) , ZMAX ) 

2 Z M I N = AM INI ( Z 2 ( I , J ) , Z M IN ) 

MAXZ=( ZMAX-ZMINJ/3.5 

IF (CUBE) GO TO 3 

MAX X = AMA X 1 (MAXX,MAXY,MAXZ ) 

MAX Y=MAXX 
MAX Z = MA XX 

3 DO 4 1=1, NX1 
DO 4 J = 1 , NY 1 

X1(I,J)=(X1( I , J)-( XMAX+XMIN)/2. )/MAXX 
Y 1 ( I , J ) = ( Y 1 ( I , J )-< YMAX+YMIN)/2. )/MAXY 

4 Z1(I,J)=(Z1( I , J)-(ZMAX+ZMIN)/2. )/MAXZ 
DO 5 1=1, NX2 

DO 5 J = 1 , NY2 

X 2 ( I , J ) = ( X2 ( I , J ) — ( XMAX + XMINJ/2. ) /MAXX 
Y2 ( I ,J ) = ( Y2 ( I , J)-( YMAX + YMIN)/2. )/MAXY 

5 Z2( I , J ) = ( Z2 ( I,J)-(ZMAX+ZMIN)/2. J/MAXZ 
RETURN 

END 
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SIJHUI MU t -'H AX I S < I , OK AHiMt- ) 
logical ok awm|- 

ni -M'S I no X{4,1 ) , Y(4, 1 ) , 7(4,1) 

?K( I , 1 'IH.D) on Tfi l 
X ( 1 ♦ l ) = 0. 

X ( 2 t 1 ) = 1.0 

X ( 3, 1) = 0. 

X { 4 , 1 ) = O. 

Y ( 1 , 1) = 0. 

Y ( 2 , 1 ) - 0 . 

Y ( 3 T 1 ) = 1.0 

Y ( 4, 1 ) = 0. 

7(1,1) = 0. 

7(2,1 ) = 0. 

7(3,1) = O. 

7(4,1) = 1.0 

K h rilKM 

1 CALL PHI (X,Y, 7,4,1 ) 

IK ( OKA M m h ) on TO 2 
KM I IRM 

2 YY = Y ( 1 , 1 ) + 4.0 
77 = 7(1,1) + 4.0 
CALL PLOT (YY,7Z,3) 

CALL PLOT ( Y(2,l )+4. , 7.(2, l)+4., 2) 

CALL SYMBOL ( Y ( 2 ,1 ) + 3 • 9 Y 1 4 3 , Z(2,1)+3.9B, .1, 1HX, 0.0, 1) 

CALL PLOT ( YY , 77 , 3 ) 

CALL PLOT ( Y ( 3 , 1 ) + 4 . , Z(3,l)+4., 2) 

CALL SYMBOL ( Y ( 3 , 1 ) + 3 . 9 7 1 4 3 , Z(3,l)+3.93, .1, 1 H Y , 0.0, 1) 

CALL PLOT ( YY, ZZ ,3 ) 

CALL PLOT ( Y ( 4, 1 )+4. , 7 (4,1) +4., 2) 

CALL S Y iY B I ! L ( Y ( 4 , 1 ) +3.97143, Z(4,l)+3.95, .1, 1.H7, 0.0, 1) 

K L T UK n 

Kivl) 


SUBROUTINE TRNMAT ( ALPHA , BETA , GAMMA ) 
COMMON /MATRIX/ TMAT(3,3) 

A=ALPHA/57. 2957795 
B=BETA/57. 2957795 
0= GAMMA/57 .2957795 
S I NA=S I N ( A ) 

S I NB=S I N ( B ) 

S I NC=S I N ( C ) 

COSA=COS( A) 

COSB=COS ( B ) 

COSC=COS(C) 

TMAT ( 1,1 )=COSC*COSB 

TMAT ( 1 ,2 )=COSC*SINB*$INA-SINC*COSA 
TMAT ( 1 , 3)=C0SC*SINB*C0$A+SINC*$INA 
TMAT ( 2 , 1 )=S INC*C0S6 

TMAT ( 2,2 ^SINC^SINB^SI N A+COSC^COS A 

TMAT (2 ,3)=S INC*SINB*COSA-COSC*SINA 

TMAT(3,1)=— SINB 

TMAT(3,2)=COSB*SINA 

TMAT ( 3,3)=COSB*COSA 

RETURN 

END 
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SUBROUTINE PHI { X , V , Z , N X t NY ) 

DIMENSION X(NX f NY), Y(NX,NY), Z(NX,NY) 
COMMON /MATRIX/ TMAT(3,3) 

C 

C 

DO 1 1 = 1, NX 
DO 1 J=1,NY 


XP = TMAT( 1 , 1 )*X ( I , J ) +TM AT { 1 , 2 ) *Y ( I , J ) + TMA T < 1 , 3 ) *Z U , J ) 
YP = TMAT {2, 1 )*X( I , J )+TMAT(2, 2)*Y ( I , J ) +TMAT ( 2 , 3 ) *Z < I , J ) 
ZP=TMAT(3,1)*X( I,J)+TMAT(3,2I*Y< I , J ) + TMA T ( 3 , 3 ) *Z ( I , J ) 


X( I , J ) =XP 
Y ( I , J ) = YP 
U I,J)=ZP 


RETURN 

END 
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APPENDIX B 


Routines for System 360 


The following are listings of the subroutines contained in the PLOT3D package for 
the IBM System 360/67. These routines are used with an on-line CDC microfilm re- 
corder. The internal routines which do the plotting are explained in the comments cards. 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE PL0T3D 
PURPOSE 

DRAW AND ROTATE THREE-DIMENSIONAL SURFACES WITH 360/67 
DD80 MICROFILM RECORDER 

USAGE 

CALL PL0T3D! XMIN,XMAX, YMIN, YMAX , S , NX PTS, NY PL S, NX PL S, NY PTS, 

F tCUBE ) 

CALL ROTATE! ALPHA, BET A, GAMMA, DR AWME) 

COMMON /LABEL/ LAB! 15) 

DESCRIPTION OF PARAMETERS 

XMIN -MINIMUM 'X* VALUE TO BE PLOTTED 

XMAX -MAXIMUM 'X* VALUE TO BE PLOTTED 

YMIN -MINIMUM *Y' VALUE TO BE PLOTTED 

YMAX -MAXIMUM •Y 1 VALUE TO BE PLOTTED 

S -ARRAY OF LENGTH 3* ( NXPTS*NY PLS+NX PL S*NY PTS ) 

NXPTS -NUMBER OF POINTS DESCRIBING LINES IN PLANES PARALLEL 
TO THE XZ PLANE 

NYPLS -NUMBER OF PLANES PARALLEL TO THE XZ PLANE 

NX PL S -NUMBER OF PLANES PARALLEL TO THE YZ PLANE 

NYPTS -NUMBER OF POINTS DESCRIBING LINES IN PLANES PARALLEL 
TO THE YZ PLANE 

F -FUNCTION SUBPROGRAM SUPPLIED BY THE USER TO DESCRIBE 

THE SURFACE TO BE DRAWN. IT MUST BE OF THE FORM 
FUNCTION F ( X * Y ) • 

CUBE -LOGICAL VARIABLE SET TO .TRUE. IF THE THREE 

DIRECTIONS ARE TO BE SCALED INDEPENDENTLY OF EACH 
OTHER* AND SET TO .FALSE. IF THE THREE DIRECTIONS 
ARE TO BE SCALED DEPENDENTLY 

ALPHA -ANGLE OF ROTATION ABOUT THE X AXIS GIVEN IN DEGREES 

BETA -ANGLE OF ROTATION ABOUT THE Y AXIS GIVEN IN DEGREES 

GAMMA -ANGLE OF ROTATION ABOUT THE Z AXIS GIVEN IN DEGREES 

DRAWME -LOGICAL VARIABLE SET TO .TRUE. IF THE ROTATED 

VIEW IS TO BE DRAWN, AND SET TO .FALSE. IF THIS 
VIEW IS NOT TO BE DRAWN. 

LAB -ARRAY OF BCD CHARACTERS USED TO LABEL THE DRAWING. 

IF LABEL IS DESIRED, SUPPLY 15 OR LESS FOUR-CHARACTER 
WORDS IN COMMON/LABEL/ LAB(15) 

REMARKS 

A CALL TO PL0T3D DEFINES AND SCALES THE SURFACE. 
SUBSEQUENT CALLS TO ROTATE ROTATE IT TO NEW FRAME OF 
REFERENCE. THE SURFACE IS NOT RETURNED TO THE ORIGINAL 
POSITION AFTER A ROTATION. THE SURFACE IS ROTATED ABOUT 
THE STATIONARY XYZ AXES IN THIS ORDER+ FIRST ABOUT THE •X' 
AXIS THEN ABOUT *Y* AND FINALLY ABOUT *Z*. 
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SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 
SCALE 
DRAW 
AXIS 
TRNMAT 
PHI 
WRITE 

A FUNCTION SUBPROGRAM SUPPLIED BY THE USER TO 
DESCRIBE THE SURFACE TO BE DRAWN* THE NAME OF THIS 
FUNCTION SUBPROGRAM MUST APPEAR IN AN EXTERNAL 
STATEMENT OF THE USER'S CALLING PROGRAM. 

IN ADDITION* PLOT3D CALLS ROUTINES DEVELOPED HERE AT 
LEWIS FOR MICROFILM PLOTTING. THEY ARE# 

LRCURV THE BASIC ROUTINE* DRAWS AN ARRAY OF (X*Y) POINTS 
AS A POINT PLOT* VECTOR PLOT* OR SPECIAL 
CHARACTER PLOT. PL0T3D USES VECTOR PLOTS. 

LRGRID SETS DP VERTICAL AND HORIZONTAL GRID LINES 

ON THE PLOTTING SURFACE. A CALL TO THIS ROUTINE 
IS PRESENTLY REQUIRED* BUT PLOT3D USES NO GRID. 
LRANGE DEFINES THE RANGES OF X AND Y IN THE USER'S 
DATA UNITS. 

LRMRGN DESCRIBES WHAT FRACTIONS OF THE PLOTTING 
SURFACE ARE TO BE USED AS MARGINS* LEFT* 

RIGHT* BOTTOM* AND TOP. 

LRLEGN PLOTS EBCDIC ARRAYS STARTING AT (X*Y) IN 
PLOTTING UNITS. 

LRCNVT CONVERTS NUMERICAL FIELDS TO EBCDIC ARRAYS 
FOR USE BY LRLEGN. 

LRCHSZ DEFINES THE SIZE OF CHARACTERS TO BE DRAWN. 

LRNON* LRNOFF CONTROL THE INTENSITY OF PLOTS. 

SUBROUTINE PLOT3D ( XMIN * XMAX * YMIN * YMAX * S *NXPTS * NYPLS * 

♦ NXPLS*NYPTS*F*CUBE) 


LOGICAL CUBE 

DIMENSION S ( 1 ) * A X ( 1 ) * BY ( 1 ) 

DATA AX * BY /O. 0*0.0/ 

FNXI = NXPTS-I 

FNY1 = NYPLS-1 

FNX2 = NXPLS-1 

FNY2 = NYPTS-1 

LI = NXPTS*NYPLS 

L2 = NXPLS*NYPTS 

N1 = LI 

N2 = Nl+Ll 

N3 = N2+LI 

N4 = N3+L2 

N5 = N4+L2 

DX = XMAX-XMIN 

DY = YMAX— YM IN 

CALL LRGRID(1*I*0. 0*0.0) 

CALL LRANGE (-3. 35* 5. 45 *-3.55 *5. 45) 
CALL LRMRGNU. 0*0. 0*3. 0*0.0) 

CALL LRCHSZ ( 2 ) 


DO 2 1=1 , NXPTS 
DO 2 J = I * NYPL S 

I I = NXPTS*! J— 1 ) + I 

II = II 

12 = Nl+II 

13 = N2+II 

S(I1) = XMIN+DX*FL0ATU-1)/FNX1 
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S ( 1 2 ) = YMIN+DY*FLOAT( J-D/FNYl 

2 S ( 13) = F ( S ( Il) t S( 12) ) 

C 

00 3 1 = ] ♦ NX PL S 
Df) 3 J=1,NYPTS • 

1 1 = K> Y P T S * ( 1-1 ) + J 

11 = N3+u 

12 = N4+II 
T 3 = NS+II 

S ( II) = XmIN + I>X*FLOAT( 1-1 )/hNX2 
S ( 12) = Y Pi I N + 0 Y * F L 0 A T { J - 1 ) / F W Y 2 

3 S ( I 3 ) = F ( S ( 11) f S ( 1 2 ) ) 

c 

c 

N 1 = N 1 + L 
w 2 = n 2 + 1 
N 3 = H3 + 1 


N4 = M4+1 ] 

i\i 5 = M 5 + 1. ] 

CALL SCALE { X m T i\i f X M A X » YM I M, YnAX . S ( 1 ) . S ( 1 ) * S ( N2 ) * n x P T S , NY PL S ♦ 

* $(N3) ,SfN4) ,S{Nb) ,NXPLS,NYPTStCUSL) 

C 1 

C 1 

CALL AXIS 10, .FALSE. ) 

C 1 

C 1 

A SUM = 0.0 I 

KSUM = 0.0 1 

CStlN =0.0 1 

CALL nki IL( ASOM, rtSUM,CS*IM) 

RETURN I 

C 1 

C I 

L iv T R Y Pill ATE ( ALPHA .BETA . G A H m A f OP A WMF ) 

LUC I CAL OR A WMF 

CALL T R NO A T ( ALPHA.PL TA. GAMMA ) 1 

CALL PHI ( S ( I ) ,S(N1 ) ,S(W2) .NXPTS .NYPLS ) 

CALL PHT ( S ( w 3 ) ,S(N4) f S(N5) . NY P T S . w X PI. S ) 

CALL AXIS (I? ORA WMF) 

A S U Pi = ASIJm+ALPHA 1 

B S 1 ) M = RSUm + KPTA 1 

CSUM = CStlN + OAHMA 1 

l r ( .NUT. DRAWEE ) PL TURN 
CALL OP A H { S { lM 1 ) ,S(M2) ,MXPTS,NYPLS) 

CALL OP A 'J ( S ( N4 ) ,S(N5> * NY PT S , NXPL S ) 

C 1 

C 

CALL HP I I PS (ASUMf HSUNfCSDM) 

CALL LRCURW(AX,HYfl ,1,1HX,1.0) 

P p I'llP N 
FNO 


BLOCK DATA 

COMMON /LABEL/ LAB ( 15 ) 
DATA LAB / 1 5*4H / 
END 
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SUBROUTINE DRAW ( Y , Z ,NPTS , NPLS ) 

DIMENSION Y ( NPT S » NPL S ) , Z ( NPT S , NPLS ) 

00 1 1=1, NPLS 

CALL LRCURVI Y ( 1 , I ) , Z(1,I> , NP TS , 2 , 1 HX , 0 . 0 ) 

1 CONTINUE 
RETURN 
END 


SUKROUT I .ME WR I TE ( Al , HI ,C1 ) 

DIMENSION AXI (3) * A X A ( 3 ) ,AYT (3) , A Y A ( 3 ) » A Z I ( 3 ) , AZ A ( 3 ) , ALPH ( 3 ) , BETA ( 3 

* ) , GAr-H 3 ) , AL IMF ( 14) ,BL I ME ( 14) ,CL I M F ( 14 ) 

DIMENSION HXI (2) , B X A ( 2 ) tBYI (2 ) ,BYA(2),HZT (2) , B Z A ( 2 ) 

COMMON /MAXES/ XMIM.XM AX, YMIM, YMAX , / MIN, /.MAX 
COMMON /LABEL/ LAR(IB) 

EOO I VALENCE ( B A , AL I NE ( 1 ) ) , ( AL PH ( 1 ) , Al. I NR ( 2 ) ) , ( KX I , AL I NE ( *> ) ) , 

* ( AXI ( 1 ) ,A|. INE<7 ) ) , (BXA, ALINE( 10) ) , ( AXA< 1 ) . ALINE( 12) ) , ( KR, 

* BL I NE ( 1 ) ) , (BETA! 1 ) ,BL INE<2 ) ) , ( RYI,BL I ME ( 3 ) ) , ( AY I ( 1 ) ,BL INE< 7) ) , 

* ( B Y A , B L I n E ( 10) ) , ( AY A ( 1 ) , BL I NE ( 1. 2 ) ) , ( BC , CL I NE ( 1 ) ) , ( G AN ( 1 ) , 

* CLINE (2 ) ) , (BZ I ,CL INEf 5 ) ) , ( AZ I ( 1 ) ,CLInE(7 ) ) , ( BZ A, CLINE ( 10) ) , 

* ( AZA(l) ,CL IME( 12 ) ) 

EOO I V AL ENCE ( EAKER,EAK IR ) 

DATA BA,BB,BC ,BX I ,BXA,BY I ,BYA ,BZ I ,BZA / 4H A = , AH B =, 4H C =, 

* 4H XM, 4HIN =, 4H XM, 4H A X =, 4H Y I'-i , 4HIN =, 4H YW, 4H A X =, 

if. 4H Z M , 4H I M =, 4H ZM, 4HAX = / 

SICK(XXX) = XXX - ELOATI I E I X ( XXX/ 340. 0 ) )*360.0 
E AKER = EAKIR 

CALL LRCMVT (XM IN, 3, AXI, 4, 12, 5) 

CALL LHCNVT (XMAX,3,AXA,4,1?,B) 

CALL LRCNVT ( YM I N , 3 , A Y I , 4 , 12 , 3 ) 

CALL LRCMVT ( YMAX, 3, AY A, 4, 12 ,5 ) 

CALL LRCNVT ( Z M I N , 3 , A Z I ,4,12,5) 

CALL LRCNVT ( ZMA X , 3 , AZ A , 4 , 1 2 , 3 ) 

RETURN 

ENTRY NRITFS ( A1,B1 ,C1 ) 

CALL L R NON 
A = SICK( Al ) 

B = SICK(Bl) 

C = SICK (Cl) 

CALL LRCMVT (A, 3, AL PH, 3, 12, 5) 

CALL LRCNVT (B, 3, BETA, 3, 12,3) 

CALL LRC NVT ( C » 3 , GAM ,3,12,3) 

CALL LRLEGN(LAB,60,0,2.0,2.66,0. ) 

CALL LRLEGNt ALINE, 36, 0,2. 0,2. 0,0. ) 

CALL LRL EGM( KLINE, 36, 0,2.0, 1 .33,0. ) 

CALL LRL EGN( CLINE, 36, 0,2. 0,0.66,0. ) 

CALL LRNOEE 

RETURN 

EMO 



SUBROUTINE SC AL E ( XM I N , XMAX , YM 1 N , YMAX , X 1 , Y 1 , Z 1 , NX 1 , NY 1 , X2 , Y2 , Z 2 , NX2 

* , NY2 *CU8E ) 

DIMENSION Xl(NXltNYl) »Y1(NX1»NY1) * Z 1 ( NX 1 » NY 1 ) , 

* X2 (NY2.NX2) , Y2 ( NY2 , NX 2 ) , Z 2 ( NY2 t NX2 ) 

COMMON /MAXES/ XM I N 1 , XMAX 1 , YM I N 1 , YMAX 1 , ZM I N , Z MAX 
LOGICAL CUBE 

REAL MAXX, MAXY, MAXZ 

XMIN1 = XM I N 

XMAX1 = XMAX 

YMIN1 = YM I N 

YMAX1 = YMAX 

MAXX = (XMAX-XMINI/3.5 

MAXY = ( YMAX-YMINI/3.5 

ZMAX = 21(1,11 

ZMIN = Zl(l,l) 

DO 1 1=1, NX1 
DO 1 J= 1 » NY 1 

ZMAX = AMAXKZ1I I, J), ZMAX) 

1 ZMIN = AMIN1 ( Z1 ( I , J ) ,ZMIN) 

DO 2 I =1 »NY2 

DO 2 J= 1 » NX2 

ZMAX = A MAXI (Z2( I , J ) , ZMAX ) 

2 ZMIN = AMIN1IZ2I I,J) , ZMIN ) 

MAXZ = (ZMAX— ZMINI/3.5 

I F ( CUBE ) GO TO 21 

MAXX = AMAX1 (MAXX, MAXY, MAXZ) 

MAXY = MAXX 
MAXZ = MAXX 
21 DO 3 1=1, NX1 
DO 3 J = 1 , N Y 1 

X1(I,J) = (Xl(I,J)-(XMAX+XMIN)/2.)/MAXX 

Y 1 ( I , J ) = (Yl(I,J)-(YMAX+YMlN)/2.)/MAXY 

3 Z 1 ( I , J ) = ( Z1 ( I , J )-( ZMAX+ZMINJ/2. ) /MAXZ 

DO 4 1=1, NY2 

00 4 J= 1 * NX2 

X2 ( I , J ) = (X2( I ,J)-(XMAX+XMIN)/2,)/MAXX 
Y2 ( I , J ) = ( Y2 ( I , J )-( YMAX+YMIN) /2. )/MAXY 

^ Z2 ( I , J ) = (Z2( I,J)-(ZMAX+ZMIN)/2.)/MAXZ 
RETURN 
END 


SIlBROUT I ME AX I S ( I , OR ANME ) 

LOGICAL URAWME 

D I nENS I ON X(4,1),Y(',,1),Z(4,1> 

DIMENSION X A ( 2 ) , X B ( 2 ) , XC ( 2 ) , Y A ( 2 ) , Y B ( 2 ) , YC ( 2 ) 
I E ( T .ME. 0) GO TO 1 
X ( 1 , 1 ) = 0 . 

X ( 2 , 1 ) = .A 

X ( 3 , 1 ) = 0 . 

X(A,1) = 0. 

Y ( 1 , 1 ) = O. 

Y ( 2 , 1 ) = 0 . 
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4 


Y ( 3 , 1 ) = 

= .4 



Y(4,l ) = 

= o. 



7(1,1) = 

= * o . 



7(2,1) = 

= o. 



7(3,1) = 

= 0. 



7(4,1) = 

= .4 



RETURN 




l CALL PHIfX,Y,J! 

, 4 , 

1 ) 

I I- ( DR A 

■IE) GO 

TO 

2 

RETURN 




2 X A ( 1 ) = 

Y ( 1 , 1 ) 

+ 

4.5 

X B ( 1 ) = 

X A ( 1 ) 



xcm = 

XA (1) 



XA(2) = 

Y( 2, 1 ) 

+ 

4.5 

X H ( 2 ) = 

Y ( 3 , 1 ) 

+ 

4.5 

XC ( 2 ) = 

Y ( 4 , 1 ) 

+ 

4.5 

Y A ( l ) = 

Z(l,l) 

+ 

3.0 

YH( 1 ) = 

Y A ( 1 ) 



Ycm = 

Y A ( 1 ) 



YA ( 2 ) = 

7(2,1) 

+ 

3.0 

Y B ( 2 ) = 

Z ( 3 , 1 ) 

+ 

3.0 

YC ( 2 ) = 

Z (4, 1 ) 

+ 

3.0 


CALL LRCURV(XA,YA,2,2,1HX,0.) 

CALL LRLArtL(lHX,l,0,XA<2) -.05,YA(2) -.05,0.) 
CALL LRCURV ( XH,YH,2,2 , 1HX,0. ) 

CALL LRLABL(1HY,1,0,XH(2) -.03,YH(2) -.05,0.) 
CALL LRCURV( XC,YC,?,2, 1HX,0. ) 

CALL LRI.AftL(]HZ,l,0,XC<2) -.Ob,YC<2) 

RETURN 
E w i ) 


SUBROUTINE TRNMAT ( ALPHA , BETA f GAMMA ) 
COMMON /MATRIX/ TMAT (3,3) 

A = ALPHA/57. 2937795 
B = BETA/57.2957795 
C = GAMMA/57.2957795 
SINA = SIN(A) 

SINB = SIN(B) 

SINC = SI NIC) 

COSA = COS(A) 

COSB = COS ( B ) 

COSC = COS(C) 

TMAT ( 1 t 1 ) = COSC-COSB 

TMAT ( 1*2 ) = COSC*SINB*SINA-SINC*COSA 
TMAT (1,3) = COSC*SINB*COSA+SINC*SINA 
TMAT (2,1) = SINC*COSB 

TMAT (2,2) = SINC^SINB^S INA+COSC^COSA 

T MA T ( 2 , 3 ) = SINC*SIN6*C0$A-C0SC*SINA 

TMAT ( 3 , 1 ) = -SINB 

TMAT ( 3 , 2 ) = COSB*SINA 

T MAT (3*3) = COSB*COSA 

RETURN 

END 
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SUBROUTINE PH I ( X , Y , Z T NX , N Y ) 

DIMENSION X ( NX, NY ) t Y< NX ,NY ) ,Z ( NX , NY ) 
COMMON /MATRIX/ TMAT{3,3) 


DO 1 I = 1 t NX 
DO 1 J = 1 f NY 


C 

c 


XP 

YP 

ZP 


T MAT ( 1 * 1 ) *X ( I f J ) + 
TMA T ( 2 » 1 ) *X ( It J) + 
T MAT ( 3 y 1 ) *X ( I t J ) + 


TMATt 1 f 2 )*Y ( I f J ) 
TM A T ( 2 f 2 ) * Y ( I , J) 
TMATI 3 t 2 )*Y ( I v J ) 


+ TM AT ( It 3)*Z<I t J ) 
+ TMAT(2f3)*Z< I f J) 
+ TMAT(3f3)*Z( I f J) 


X( If J) = XP 
Y( I , J ) = YP 
1 Z ( I f J) = ZP 


C 

c 


RETURN 

END 
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Figure 1. - Function from statistical mechanics, z - e"^ x+ *^*(xyf 3,z . 
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Figure 2. - Surface fit to points derived from experiment and analysis combined. 
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Figure 3. - Particular solution to the two-dimensional wave equation, 
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Figure 4. - Function with two poles of different order. 
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(b) View 2. 
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Figure 4. - Continued. 
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(c) View 3. 
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Figure 4. - Concluded. 
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Fixed axes Figure axes (A is rotation matrix for a, see ref. 10) 

(a) Rotation about x-axis. 
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(b) Rotation about y-axis. 
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(c) Rotation about z-axis. 

Figure 5. - Rotations of the surface. 
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Figure 6. - General flow chart for PLOT3D. 
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